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We propose that electron doped nontransition metal-phthalocyanines (MPc) like ZnPc and MgPc, 
similar to those very recently reported, should constitute novel strongly correlated metals. Due 
to orbital degeneracy, Jahn- Teller coupling and Hund's rule exchange, and with a large on-site 
Coulomb repulsion, these molecular conductors should display, particularly near half filling at two 
electrons/molecule, very unconventional properties, including Mott insulators, strongly correlated 
superconductivity, and other intriguing phases. 
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Building novel metals by doping molecular crys- 
tals such as polyacetylene, fullerenes, TTF-TCNQ, 
(TMTSF) 2 X, (TMTTF) 2 X and (BEDT-TTF) 2 X salts, 
etc., is a well trodden route [l|,|2j, but remains an excit- 
ing and ever moving front. Very recently the Delft group 
showed that thin films of transition metal phthalocya- 
nines (MPcs) FePc, CoPc, NiPc, CuPc, initially insulat- 
ing, can be turned genuinely metallic through potassium 
doping Q ■ Electron doping appears to take place largely 
in the twofold degenerate lowest unoccupied e g molecu- 
lar orbital (LUMO) [111] of the molecules. The simplest 
rigid band model naturally explains, according to [3j , why 
the MPcs, initially insulating when pristine (n = 0) [6J], 
become metallic upon increasing doping (0 < n < 4), 
ending up again as insulators at full doping (n = 4). It 
is not yet clear whether stoichiometric compound phases 
may exist here as they do in alkali doped fullerides, but 
there is an otherwise striking overall analogy, to the point 
that even the conductance values reported at optimal 
metallic doping are close in magnitude and temperature 
(in)dependence to those of K^Ceo films. 

While the pursuit of metallicity-cum-magnetism in the 
transition metal doped MPcs |5j will in the future con- 
stitute an interesting goal in itself, the scope of this note 
is to point out newer directions and possibilities that 
can make slightly different metallic doped MPcs, yet 
to be realized, potentially even more exciting. Doped 
MPcs arc not, we claim, regular metals, but consti- 
tute strongly correlated electron systems, akin to doped 
two-band Mott insulators. Characterized by a strong 
on-site electron repulsion, by the e g orbital degener- 
acy, and by intra-site Jahn- Teller and electron-electron 
multiplet couplings, the doped MPc molecular system 
should approach near n = 2, as foreshadowed by re- 
cent studies 0, la Oj a novel unstable fixed point herald- 
ing a wealth of possible low-temperature phenomena and 
phases. We conduct in this letter a preliminary explo- 
ration of this scenario by addressing theoretically - and 
thus proposing the experimental realization of - new met- 



als obtained through electron doping of ZnPc, MgPc 
and other such non-transition metal MPcs. The extra 
electrons should flow into the MPc 2e g lowest unoccu- 
pied molecular orbital (LUMO)|4|, to form a two band 
metal. The LUMO states are ligand shell orbitals loosely 
surrounding the central metal ion, with a large intra- 
molecular Coulomb repulsion U in comparison with the 
narrow electron bandwidth W expected from weak inter- 
molecular electron hopping. In the neighborhood of half 
filling, n = 2, a doubly degenerate e g orbital also has a 
molecular Hund's rule exchange J (lowering the triplet 
state by 4|J| relative to the singlet), and a Jahn- Teller 
(JT) coupling of the electronic state to B lg and B 2g MPc 
molecular vibrations, lowering the singlet state energy by 
an amount Ejt- Thus there will be some cancellation be- 
tween the two couplings. These ingredients are similar in 
nature to those present in fullerides, where orders of mag- 
nitude are U = 1 eV, E JT = 0.15 eV, J = 0.03 eV, and 
IF = 0.5 eV. Although the corresponding values in MPcs 
are not yet so accurately known, we obtain estimates that 
are surprisingly close. The similarity noted in Ref. |3j be- 
tween alkali doped MPcs and fullerides |2| may thus reach 
much deeper than just the dopability and the closeness 
of mobilities pointed out experimentally. As in fullerides, 
the resulting phase diagram will depend upon the bal- 
ance between low and high spin states of the molecular 
ion. Unlike fullerides, however, the twofold degeneracy 
will reveal a potentially richer overall phase diagram. 

The Hamiltonian for electrons in the e g -orbitals of 
molecule i is 
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where hi is the overall occupation at site z, fi ia , with 
a = x,y,z, are pseudo-spin operators spanning the 
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where r a are Pauli matrices. Here c\ ary denotes creation 
of electrons with spin a at site i in orbital a = 1,2 of 
the LUMO e g doublet. The last two terms in (|TJ de- 
scribe JT coupling of strength g to a doubly degenerate 
vibration with frequency TUoq, coordinates q^i and 5^2 
and momenta pi t i and pi^, where Ejt — &g 2 /fru>o- A 
single but degenerate mode is assumed as a simplifying 
approximation, replacing the actual 14 pairs of nondc- 
generate JT modes B lg and B 2g allowed by the MPc 
D4/1 symmetry. As detailed in note ll(J| an estimate of 
parameters appropriate to MPcs such as ZnPc, MgPc in- 
dicates U ~ leV and, for n = 2, a singlet energy gain 
Ejt -4| J| - 0.06-0.07 eV over Hund's rule triplet. This 
prevalence of singlet over triplet, the main ingredient po- 
tentially leading to s-wave superconductivity, could be 
mitigated or even reversed in a solid state compound; we 
shall therefore examine both possibilities in the rest of 
this paper. 

In the MPc molecular crystal, electrons can quantum 
mechanically hop between neighboring sites 
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Assuming first a generic three-dimensional structure with 
v neighboring molecules and a weak isotropic first- 
neighbor hopping t, the LUMO e g orbitals will give rise 
to a pair of narrow molecular bands of width W = 2\t\v. 
(We will discuss the additional signature of a quasi-one- 
dimensional structure of these compounds in the final 
part of the paper.) Before doping we generally expect a 
narrow band W ~ 0.3 eV (see below), due to very weak 
van der Waals bonding between molecules. While W 
will eventually depend on the intermolecular geometry, 
including the way that may be altered by doping, we will 
assume W to remain of approximately the same magni- 
tude after doping. In addition to (JTJ and (J2J wc should in 
principle include additional terms such as the dispersion 
of vibrations (to account for the coupling between the JT 
distortion of a molecule and those of other molecules), 
inter-molecular electron-electron interactions, and more. 
We will presently neglect them, concentrating on the ba- 
sic Hamiltonian 
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Because U/W is large, and for n ~ 2, this system is gen- 
erally close to a Mott transition. There, the quasiparticlc 
band is very narrow and it is a good approximation to 
neglect retardation, replacing (JJJ with 
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where J e // = J — Ejt/4. 

The further simplifying assumption of nearest- 



neighbor hopping that are diagonal in the orbital in- 
dices |ll| brings this Hamiltonian to coincide with 
that rece ntly studied by Dynamical Mean Field Theory 
(DMFT)[jj in Refs. 00. Close to n = 2 it was shown 
to display a rich phase diagram including, besides the 
regular metal, a Mott insulator, a strongly correlated su- 
perconductor, and a pseudo-gap metal. Within DMFT, 
the lattice model (0 is mapped onto an Anderson im- 
purity self-consistently coupled to a conduction electron 
bath [l2j . The many-body physics of the Anderson im- 
purity and its fixed points foreshadow that of the infi- 
nite many-body system[9| , as confirmed by the calculated 
DMFT phase diagram [7), 18J . We exploit here these re- 
sults to analyse the implications for the phase diagram 
of half-filling doped MPcs considering both possibilities, 
namely: (i) singlet ground state or (ii) triplet ground 
state. 

(i) J e ff < - The n—2 molecular ground state is a 
non degenerate spin singlet accompanied by a dynamical 
Jahn- Teller effect. If in the solid U/W is larger than the 
critical Mott value (between 1 and 2, depending in detail 
on J and Ejt) the half filled MPc will realize a non- 
magnetic, singlet Mott insulator |l3j. For U just below 
that value, or for light doping away from n=2, the MPc 
solid will be lightly metallic. In this "doped Mott insula- 
tor" regime the Anderson impurity displays instabilities 
against symmetry-broken phases in both particle-particle 
and particle-hole channels; which instability will even- 
tually win out depends on band-structure and coupling 
details. If nesting or other band structure singularities 
are absent or very weak, the particle-particle instabil- 
ity should dominate leading to an s-wave superconductor 
with order parameter 
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This state is a strongly correlated superconductor (SCS) 
of a kind first pointed out in Ref. [bj|, further con- 
firmed in the two band model in Refs. [7J, |8|- It was 
also suggested |7J, [lj] - although without a calculation 
of T c - that SCS is a "high temperature superconduc- 
tor" in the sense that the SCS superconducting gap may 
reach values several orders of magnitudes greater than 
the corresponding BCS value (that would be attained for 
U = 0), A BC s ~ hui e~ 1/x , where A = 2p \J ef f\, p 
is the density of states at the Fermi level per spin and 
band. 

If, on the contrary, the particle-hole instability channel 
was favored by some detail of the band structure such as 
Fermi surface nesting, then one of two alternative symme- 
try broken phases can be expected. The first is a trivial 
one, with order parameter 

Aj T = cos^(f x )+sin(?!>(fj / ) ^0, (6) 

corresponding to a cooperative Jahn- Teller distorted 



state, either modulated or uniform, the angle cf> reflect- 
ing U(l) orbital symmetry breaking. The model in fact 
possesses 0(2) orbital symmetry, which naturally de- 
composes into Z 2 ® U(l). Z 2 is the discrete symmetry 
t z — > — t z , while U(l) represents invariance under ro- 
tations around the pseudo-spin z-axis. As a consequence 
of the static cooperative Jahn- Teller distortion, a con- 
ventional insulating band gap or at least a strong low- 
ering of metallic density of states should take place at 
the Fermi level. More interesting is the alternative insta- 
bility, associated with a different order parameter in the 
particle-hole channel 

Apt = ^2{c\ a ct q/ 3 ■ Sc 1/3 ) - (c| Q <7 Q/3 ■ S c 2f} ) ^ , (7) 

a/3 

breaking both orbital Z 2 and spin SU(2) symmetry, S 
being the direction along which the spin SU(2) symmetry 
is broken. We note that the LUMO e g orbitals are mainly 
localized onto the MPc four N-atom ring and have odd 
parity with respect to the central metal ion. Therefore 
the discrete Z 2 symmetry translates into parity in MPcs, 
whence Apt ^ represents a spin current flowing on 
the Pc ring, a very intriguing possibility. This kind of 
phase should be either insulating or a poor metal like 
the cooperative JT phase. Unlike the latter, it should 
break both time reversal and parity, conserving only the 
product of the two. These properties are reminiscent of 
a phase discussed some time ago by Varma [15j . 

(ii) J e ff > - In this case (less likely in MPcs than (i), 
as said above) the isolated n =2 molecular ion is a spin 
triplet. The strongest instability here is towards bulk 
magnetism 9], ferro or antiferro depending on the band 
structure, a state which may or may not be accompanied 
by parity-symmetry breaking or by a Jahn- Teller distor- 
tion. If (depending on lattice structure) magnetism was 
sufficiently frustrated, then it would be possible for spin- 
triplet superconductivity to appear with order parameter 
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Since orbitals 1 and 2 have opposite parity, this translates 
into a p-wave spin-triplet superconducting order param- 
eter. This is an interesting possibility since, at variance 
with conventional p-wave superconductors, the pairing 
function is nonzero on the molecule. According to the 
DMFT analysis of Ref. 8], where this kind of triplet 
suoerconductor was first discussed, this superconducting 
instability is also enhanced by the proximity of a Mott 
insulating phase, although not as dramatically as in case 

CO- 

Much of the uncertainty in the above scenario depends 
on our ignorance of the doped MPc film crystal structure, 
so far treated generically. A second guess based on the 
crystal structure of pristine MPcs suggests that the al- 
kali doped crystals could, similar to the undoped alpha- 




FIG. 1: Density functional electronic band structure of un- 
doped a structure MgPc. The molecular stack direction is 
r — Y. The crystal structure was taken from Ref. |20|. The 
monoclinic cell (a = 26.29, b = 3.818, c = 23.92A, j3 = 94.6) 
contains four molecules and has symmetry C2/n. 



phase of the CuPc films |3(], be made up of quasi-one- 
dimensional (ID) chains. We carried out density func- 
tional electronic structure calculations of MgPc assuming 
an alpha phase structure with n = 16], and found an 
intra-chain e g bandwidth W ~ 0.3 eV but an inter-chain 
bandwidth between one and two orders of magnitude nar- 
rower (see Fig.p). If doped MPcs were indeed this close to 
ID, then our earlier DMFT-based analysis, valid for the 
opposite limit of large coordination lattices, should be re- 
placed by another where the additional effect of quantum 
fluctuations and other ID specific anomalies are properly 
treated. 

We studied this ID model by adapting the two-loop 
renormalization group (RG) equations derived by Ref. 
|17| to the two-band model. The fixed point Hamilto- 
nian towards which the model flows under RG was an- 
alyzed by means of bosonization|l7J, |l8|. At strict half- 
filling we found a spin-gapped insulator for any value of 
U. However, for large U, this was a spin-liquid Mott in- 
sulator, i.e. a spin-gapped phase without any symmetry 
breaking. That corresponds either to the DMFT singlet 
Mott state found with J e ff < 0, or to a Haldane spin-1 
chain for J e ff > 0. Here, the DMFT metallic phase at 
small U is replaced by a spontaneously dimerized insu- 
lator, i.e. a spin-gapped phase with broken translational 
symmetry, for either sign of J e ff ■ Away from half-filling 
the model turns metallic: nonetheless we find that the 
spin-gap survives. In particular, in the large U regime 
likely pertinent to the MPcs, the spin-liquid Mott insu- 
lator at J e ff > turns at low doping into a spin-gapped 
metal with dominant fluctuations in the SC singlet chan- 
nel (J5J and sub-dominant fluctuations in the particle-hole 



Akp channels identified by Ref. Q- The weakly-doped 
spin-1 chain at J e ff < feels on the contrary more 
the effects of the reduced dimensionality which causes 
the spin gap to survive in the metallic phase, hence pre- 
venting long-range spin-triplet superconducting correla- 
tions. Indeed, in ID the leading SC fluctuations appear 
in a spin-singlet orbital-singlet but space-odd particle- 
particle channel |l9|. 

Weak but nonzero interchain interactions will in a hy- 
pothetical doped MPc with alpha phase-like structure 
turn all these instabilities into properly long-range or- 
dered phases at sufficiently low temperatures. While our 
present state of ignorance prevents further in depth dis- 
cussion, the expected emerging scenario is parallel to that 
obtained by DMFT, apart from those differences specific 
to ID, like dimerization at n — 2 or the persistence of a 
spin gap away from half-filling for regular Hund's rules. 

Summarizing, we conclude first of all that stoichiomet- 
rically doped nontransition MPcs should display Mott in- 
sulating phases at all integer fillings, in particular at half 
filling. Analysis of a simple model 0J revealed a stun- 
ning variety of phases in the immediate neighborhood of 
the Mott metal-insulator transition near half filling. The 
variety is greater than either that predicted for the three 
band case fullerides 14], or that experimentally known 
in fullerides |2j as well as in doped organics (with lower 
symmetry than MPcs) near half filling[l|. 

Many of the phases described above will individually 
merit an in-depth study. At the theoretical level, ap- 
proximations will need to be improved, for example by 
treating properly the non-degeneracy of the Bi g and _B 2ff 
molecular vibrations, the non-diagonal hopping, the dis- 
persion of vibrations, etc. At the moment, nonetheless, 
it seems fair to say that more direct experimental input 
is urgent and essential. It will be important to dope 
and metallize non-transition MPcs, and to investigate 
whether they can be made at least locally stoichiomet- 
ric. Subsequent examination by, e.g., scanning tunneling 
spectroscopy could search for the Mott insulating state 
and begin broaching in its neighborhood this fantastic 
scenario. 
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